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Abstract 

This paper is devoted to the study of the regularity of solutions to some systems of reaction- 
diffusion equations, with reaction terms having a subquadratic growth. We show the global 
boundedness and regularity of solutions, without smallness assumptions, in any dimension 
N. The proof is based on blow-up techniques. The natural entropy of the system plays a 
crucial role in the analysis. It allows us to use of De Giorgi type methods introduced for 
elliptic regularity with rough coefficients. In spite these systems are entropy supercritical, it 
is possible to control the hypothetical blow-ups, in the critical scaling, via a very weak norm. 
Analogies with the Navier-Stokes equation are briefly discussed in the introduction. 
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1 Introduction 

This paper is dedicated to the study of the global regularity in time of the solutions to a class of 
react ion- diffusion systems. Reaction-diffusion systems are used as models for a variety of problems, 
especially in chemistry and biology [6j HU [131 El EH EH] . The question of the existence of global 
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solutions is particularly important and it has been widely studied [201 E21 I2U ED]- In full 
generality, such system can have solutions which blow up in finite time as it is well known 
when considering non linear heat equations [T5l 151] . In this article we focus on systems which 
are entropy supercritical, that is, such that the preserved physical quantities, as the mass and the 
entropy, are not shrinking (not even preserved) via the universal scaling of the system. The global 
regularity for entropy (or energy) supercritical problems is an important question in several areas 
like nonlinear waves, Schroedinger, and the Navier-Stokes equations in dimension N > 3. 

We present, here, a class of systems for which we are able to show global regularity despite their 
supercritical nature. 

We consider systems of the form: 

d t a - DAa = Q(a), t > 0, x G R N , 

(1) 

a(0,x) = (o?(x),a§(z),- • ;a° P (x)), x G R , P > 1. 
We assume that the matrix of diffusion is diagonal: 

Di G R, for 1 < z < P, 
and d := sup {Di}. 

l<i<P 

In addition, we assume that the reaction diffusion term Q is regular and it satisfies the following 
four conditions: there exists < v < 2 and A > such that 



D = diag(D 1 ,...,D P ), 

with 

d := inf {DA > 

— Ki<P 



Qi(a) > 0, if ai < 0, 1 < i < P, (2) 

\VQi(a)\ < A\a\ u -\ aGM+, 1 < i < P, (3) 
p 

^Q,(a)=0, a G (4) 
p 

^ In a<Qi(a) < 0, a G Rj. (5) 



i=l 



Hypothesis (J2J) ensures the nonnegativity of the a,, hypothesis ([5]) is a restriction on the growth 
of Q, hypothesis (jlj) ensures the conservation of the total mass Yld=i I a i & x -> an d hypothesis ([S]) 
ensures the non increase of the entropy Y2i J a i ^ n a i dx. Note that (J3j is the only restrictive hy- 
pothesis. It requires that the system is at most subquadratic. The other hypothesis are physical 
and pretty standard in systems coming from chemistry. 
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Examples of systems verifying such hypothesis for P = 4 are given by: 

Qi(a) = (-l)*(0(aia 3 ) - 0(a 2 a 4 )), a = (ai, a 2 , a 3 , a 4 ) G R+, 1 < i < 4, 
where G C°°(R + ) is any nondecreasing function verifying 

(f)(z) = 0, for z < 0, 
0(z) = ^/ 2 , for Z >1. 

We consider nonneg ative initial values a° G C°°(R N ) n L 00 ^), 1 < i < P, verifying 

sup < / a°(l + |sc| + |ln(a-)|) dx > < oo, (6) 

but without smallness condition. We note that these additional constraints correspond to the finite 
initial mass and total entropy, together with a condition of confinement of the mass near the origin. 

The main result of this paper is the following theorem. 

Theorem 1 Let u < 2. Consider a system (T7]j verifying hypothesis |D) ; (T3j) ; Q) ; and (Tjp. Then, 
for any a = (a?, - ■ -,a° P ), with a° G C°°(R N ) H 7°°^), a? > /or 1 < i < P, verifying the 
condition i/iere exists a unique smooth global solution a = (ai, • • •, ap) defined on [0, oo) x R . 
Moreover, > /or a// 1 < z < P. 

A system of particular interest, which is a limit case (not included) of this study, is the quadratic 
system of reaction with four species corresponding to the chemical reaction 

At + A 3 ^ A 2 + A A . 

In this case the system (JT|) has P = 4, and 

Qi(a) = (-l) 4 (aia 3 - a 2 a 4 ), a = (ai, a 2 , a 3 , a 4 ) G R+, 1 < z < 4, 

where a« is the concentration of the species A4. Let us denote 

Ui = y/ai, 1 < i < 4. 

The mass conservation and entropy dissipation give (see section [2]) the following bounds for any 
T > 0: 

17 G 7°°(0, T; 7 2 (lR Ar )), V S 17 G 7 2 ((0, T) x R^). 
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Also, the universal scaling is given by 

U e (t,x) = ell(e 2 t,ex) 1 e > 0. 

This means that U is a solution to the system ([T]) if and only if U £ is a solution to the same problem. 
Note that U has exactly the same conserved quantities and the exact same universal scaling as the 
Navier-Stokes equation, which is supercritical for N > 2. We interpret this by saying that our result 
gives the full regularity for a family of parabolic systems, in any dimension, which has almost the 
same supercriticality as the Navier-Stokes equation. 

However, we observe that these systems of reaction-diffusion have some features quite different 
from fluid mechanics. In the inviscid case where = for all z", 1 < i < P, a maximum principle 
ensures that the solutions remain globally bounded. But this structure breaks down in the diffusive 
case provided that the diffusion coefficients D^s are not all the same. Surprisingly, such systems 
are, actually, destabilized by the diffusion. 

One of the key arguments in our proof is the use of a local parabolic regularization effect. This 
argument basically claims that if the nonlinear terms are somehow controlled and small in a cylinder 
(to — I, to) x B(xo, 1), then the value of the solution at (to,xo) is bounded by 1. This argument 
is at the heart of the partial regularity results for the Navier-Stokes equations jU [T71 [28j [29] . 
In [31], a new proof is proposed based on the De Giorgi's parabolic regularity method [7]. This 
technique is particularly powerful and gives important results in different physical areas, such as 
the quasigeostrophic equation [5], and the compressible Fourier-Navier- Stokes system [18j[T9]. This 
method has been used for the first time in the context of reaction-diffusion equation in [TB]. This 
paper contains the cases N <2, even with quadratic growth. It corresponds to the entropy critical 
case. 

The local parabolic regularization results are based on recursive controls, along a family of 
shrinking cylinders, of the nonlinear terms in a cylinder by the entropy (or energy) on a bigger 
one. We note that, in the context of react ion- diffusion systems, the classical methods based on 
the Green function of the heat equation do not work in the supercritical cases. The main problem 
is that we have different diffusion coefficients. In this context, the De Giorgi method proves to 
be particularly powerful as it exploits the physical entropy quantity jSJ). Using this property, this 
method depletes the nonlinearity of one exponent. This is one of the key facts which allows us to 
work with subquadratic nonlinearity for any dimension N. It follows that, locally, it is enough to 
control recursively slightly bigger (to the log) norms that the L 1 norm to get the regularization 
since such norms are controlled by the entropy. 

To get the global regularity, we shrink these local norms, about any point {to, x ) through the 
universal scaling 

a £ (s, y) = e~^a{e 2 s + t , ey + x ), for e > 0. 

If a = (ai, • • •, dp) is a solution to the system, then a £ = (a\ , • • •, ap) is a solution to a system with a 
reaction term Q £ which shares the same properties as the reaction term Q [Hypothesis (j2J), (EJ), PJ, 
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and (jSJ)] with the same v and the same constant A. If there is an e > 0, small enough, for which the 
norm is very small, then we can use the local parabolic regular izat ion result on a 6 = (af, • • -,ctp) 
to ensure the regularity of a = (aj, • • •, ap) at the point (to, xq). In the entropy supercritical cases, 
all the quantities obtained through the entropy blow up along the rescaling. But, on the other 
hand, we can show than the solution is globally bounded in a very weak space which has the same 
homogeneity as L oc (W~ 2,oc ). The norm in this weak space shrinks in the subquadratic growth 
v < 2 as 

£ II 4-2i/ 

a IIl 00 ^ -2 ' 00 ) = £ II oIIl 00 ^- 2 ' 00 )- 

It is the nonnegativity of the solutions that allows us to control the depleted nonlinear term by this 
weak norm. We can use, then, the local parabolic regularization property. 

The control of the solution in this weak space is based on a nice duality argument introduced 
by Pierre and Schmitt [26j. In |26j, they obtain a global bound in L 2 (L 2 ) of the solution. A similar 
argument has been used in [TD] to get global weak solutions. The L 2 (L 2 ) norm does not shrink 
through the universal scaling for the entropy supercritical cases. However, we will show that the 
same method can be applied for the weak norm that we present. 

The study of this weak norm was suggested by the paper [32], where the L 00 (W /_1 ' 00 ) norm is 
shown to play an important role in the Navier-Stokes equation. We observe that the L°°(W~ 1,0 °) 
norm has the same scaling as the ^(BMO" 1 ) norm studied by Koch and Tataru for the Navier- 
Stokes equation. The weak norm that we introduce corresponds, formally, for the quadratic case, 
to the L^iW' 2 ' 00 ) norm. 

The additional control of the weak norm comes from a linear equation verified by the total mass 
p = Y^i a i- Pierre and Schmitt [26] use this extra equation to get the L 2 (L 2 ) estimates. They also 
show that the full regularity cannot come solely from this argument as they provide some explicit 
examples of solutions to this type of linear equation which blow up in finite time. 

We note that, in the contest of the Navier-Stokes equation, an equation on the vorticity is also 
provided as an extra. For the Euler equation (the inviscid case), the vorticity u is a solution to the 
equation: 

d t uj + (u ■ V)u - (u ■ V)w = 0. 

If we drop the dependence of vorticity on the velocity u, the above equation can be seen as a linear 
equation on to, which can be solved in the Lagrangian coordinates as 

u{t,a) = (uj(0,a)-V a )X{t,a), 

where X is the flow given by 

d t X(t,a) = u(t,X(t,a)). 

Note that, for any smooth initial values, this gives a uniform control of the vorticity in L 00 (H // 1 ~ C 1 ' 00 ) 
in the Lagrangian coordinates. This norm is, actually, shrinking through the universal scaling of the 
Navier-Stokes equation. Unfortunately, this structure seems to be destroyed by the viscosity term. 
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Such proprerty provides another analogy with the reaction-diffusion systems: as for the maximum 
principle for the reaction-diffusion systems, a crucial supercritical structure known on the Euler 
system is destabilized by the diffusion term. 

Global existence of weak solutions of ([!]), was established in [10]. The dissipation property ([5]) 
is also the basic tool for studying the asymptotic trend to equilibrium [8j [9] in the spirit of the 
entropy/entropy dissipation techniques which are presented e.g. in [33] (we refer also to [2] for 
further investigation of the large time behavior of nonlinear evolution systems using the entropy 
dissipation). Let us also mention that ([T]) can be derived through hydrodynamic scaling from 
kinetic models, see [3]. 

As remarked above, the results of Theorem [1] are trivial if the diffusion coefficients TVs are all 
equal. Another trivial case corresponds to P = 2 (two species) where a maximum principle holds 
(even without the subquadratic property ([3])). For the sake of completeness we will give a proof of 
this result in the appendix as we did not find this case in the literature. 

Theorem 2 Consider a system with P = 2, and verifying hypothesis |I)) Q), and |3]). Then, 
for any a = (a?,a°) ; with a\,al G C 00 ^) nL°°(R N ), a?,aij > 0, verifying the condition there 
exists a unique smooth global solution a = (01,02) defined on [0, 00) x M. N . Moreover, a 1; a2 > 0, 
and 

sup {ai{t, x), ci2(t, x)} < sup {ai(x), a^(x)}. 

(t,x')e(0,oo)xM JV xeR N 

The proof of the maximum principle collapses for P > 3. 

Some of the estimates in this paper will be based on results obtained in [IB] , but, for the sake 
of completeness this paper is self contained. 

For any smooth initial data, standard theory ensures the existence of a smooth solution on (at 
least) a short time [0,T), T > 0. We set T to be the biggest of such lapse of time. Our aim is to 
show that T = 00. Standard bootstrapping arguments give that if T < 00 then 

lim sup ||ai(t, -)||l°° (r jv) = +00. 

t^T i<i<p 

We will obtain a uniform bound on [0, To) contradicting the blow-up of the solution in finite time. 
Indeed, we will show that this bound depends only on T an d the quantity M defined by 

M = sup I / a i {l + \x\ + \\n{a i )\)dx+\\a i \\ L o O{R N ) \ . (7) 
l<i<P Urn J 
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In the second section, we provide standard a priori estimates based on the mass conservation 
and the entropy dissipation. The third section is dedicated to the local parabolic regularization 
principle. The duality arguments are given in the fourth section. We introduce and then apply the 
rescaling arguments in the last section. 



2 Entropy Dissipation 



In this section, we derive a priori estimates on [0, T ) where the solution is smooth. The dimensional 
cases of N < 2 have been already studied in [16] , but they can also be deduced from the techniques 
that we present in this paper. (The reader can easily replace the classical inequalities in the sequel 
of the proofs for dimensions N < 2). Henceforth, we will analyze here only the dimensions N > 3. 

We discuss the a priori estimates that can be naturally deduced from (jlj) and (jSJ). 

Proposition 3 There exist two constants Co and C\, such that the following is true. Let (T7|) be 
any system verifying (TJj), @) ; and and any initial values a° > ; smooth, verifying 

M = sup (||a?|| L co (R * ) + / a°(l + \x\ + |ln(a?)|)dx J < oo. (8) 

l<i<P V JR N J 

Let a = (a 1; • • -,ap) be the associated solution. Assume that a is regular on its maximal lapse of 
time [0, T ). Let 

p 

Then, for any < T < T , 

P ~ P „T 

c * 

0<t<T 



sup I V / aj(l + \x\ + \ ln(aj)|)(t, x) dx \ + 2c? / / \V x^/aTif dx ds 
<t<T < ~^ Jrk > ~[ Jo Jr n 

P t 

+ V / / ®{a)dxds < {Co + C{T){Mo + 1). 
~l Jo Jr n 



Proof. First, thanks to (|2J), we have a^t,^) > for any t e [0,T ), a; G R^, and 1 < i < P. As a 
consequence of (jlj) and (jSJ), we get 



(it 

i=i 



2 / a>i dx = 



and 

d_ 
~dt 



Also, 



V/ ai(l + ln(ai))dx ) +Y] [ AV x a { ■ cb + / Ddx = 0. 

^"^^ y "^ii^ a i Jr n 

/ DiV x ai ■ —!— - dx > d [ — — — dx = 4d / I V xV /oi | 2 
^ ii N a» ~^ Jr n a i ~( Jr n 



Moreover 

p 



d 
eft 



/ a,i\x\ dx = — / DiV x ai-^— dx 

~^JR N ~l Jr n \ x \ 

^ f ^ f I ^ ^ I 

< rfN / |V x aj| dx = d} j / — ^=-^/aidx 

~l Jr n ~[ Jr n V®i 

— o ( dx + ^— / ai dx. 

i=l a ' 1 — i=l 



Let < t < T. We integrate on the interval [0, t], 

p r a p rt 



> / Oj(l + |x| + ln(aj)) dx + = ^, / / — - - dxds + / / Ddxds 

~l Jr n 2 ^ y y R iv aj y y R iv 

< M + — ^ j didx ds 



i=1 Jo JR N 

< (l + td 2 / (2d)) M . 
Finally, we estimate the negative part of the ajln(aj): 



/ ai\\n(ai)\dx = / a { ln(aj) dx - 2 / a i ln(a i )(l < o .< e -|x|/2 + l e -H/2<oi<i) ^ 

JR^ JR^ JR/v 

x|oj dx 

Jrjv e y M jv y M jv 

since — sln(s) < \\[s for any < s < 1. Combining together all the pieces we conclude the 
statement of the proposition. 
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3 Local parabolic regularization principle 

We consider solutions of ([1]) that are defined for negative times. The following proposition is the 
main result of this section. 

Proposition 4 Let (Q]j be a system verifying (U|), (Tjj) ; Q) ; and ([3]). Let a = (ai, ■ ■ -,ap) be a 
solution to (CJj on (—3, 0) x -8(0, 3), such that are nonnegative functions for 1 < i < P. Then, 
for any p > 1, there exists a universal constant 5* > (depending on the constant A as in (T3|) ; d, 
d, P, N and p), such that, if a = (oi, • ■ -,op) verifies 



then, < a;(0, 0) < 1 for 1 < i < P. 

The value at (0, 0) can be hence controlled by any L p norm, p > 1, on a surrounding cylinder. Such 
results is quite surprising, since we work with subquadratic reaction terms. 

In the spirit of the Stampacchia cut-off method, L°° bounds of solutions of certain PDEs can be 
deduced from the behavior of suitable non linear functionals. Here, such functionals are constructed 
in a way that they use the dissipation property ([5]). This is based on the De Giorgi techniques and 
it is reminiscent of the method introduced by Alikakos pp. 

Let us consider the non negative, C 1 , and convex function 



p 





if z > 0, 
if z < 0. 



We study the evolution of the "entropy at level K 




for R > 0, where the set B n will be specified later. 



Lemma 5 Let 1 < v < 2. For any a G M p ; and for any < R < 1, we have 

p 



\Qi( a ) - Qd 1 + l a - R }+)\ < 2PA I X + l a - R ] 



i=l 



where A is given by (TJ|) and 1 + [a — R] + := (1 + [a\ — R] + , ••-,!+ [dp — R)+)- 
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Proof. For \u\ < \v \, by the sub-quadratic growth for the reaction term Q, hypothesis ([3]), we have 

p 

~ Q<( v )\ ^ Ap \ v ~ "IM"" 1 - 

i=l 

The inequality follows by substituting (u, v) by (a, 1 + [a — R]+). ■ 

Let us set k n = 1 - l/2 n , t n = 1 + l/2 n , £„, = 5(0, t n ), Q n = 0) x B„. Note that B n C £ n _i 
and Q n C Q n _i. 

We introduce the cut-off functions 

C n : M N -> R, < Cn(a?) < 1, 

C„(x) = 1 for x G <S n , Cn(a^) = for a; G Cfi n _i, 

sup |9?-£ n (x)| < C 2 2n with C universal constant. 

i,je{i,...,N}, xm N 

Next, we give an estimate on the local dissipation of entropy at the level R. 

Proposition 6 There exists a universal constant C (depending on d,d, A and P), such that for 
every a = (a 1; • • •, ap) solution of (QP, for any < R < 1, we have 

p p 
sup {Y^ / $(aj - R)(t, x) dx] + d,y~ K j / / + [a* - i?] + | 2 dxdT 

p 

2 2n(j I I (1 + [oj - R} + ) ln(l + [a* - R]+) dxdr. 



< 



Proof. We multiply (JTJ by C n $'(ai ~ an d we sum 

d P f 

j t J2 ^ $ ( a * ~ ^) dx = A + s 

1=1 "'^n-l 



(9) 



where 



and 



p r 

i=l ■'Bn-l 

B -=Y, Q l (a)&(a l -R)( n dx. 

i=l J Bn-l 
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We rewrite A as, 
where 



A = —E + F 



and 

Moreover, 

p 



p r 

E :=J2 A|Va,| 2 - R)( n dx 

i = l JBn-l 

p r 

F := V / Di$( ai - R)A( n dx. 
<=i 



£1 ^ y^-i i + k - 



V / AV,(1 + [a, - ■ V x (l + [a, - + ) — — 

i=1 ^Sn-l 1 + i' 



>k 

dx 



di — R\. 



* f \V x (l + [g t -R] + )\ 2 



i=l JB "-i 

p 



> 



l + [ai- R]. 
4^y^ / |V x a/1 + [a* - R] + \ 2 dx. 



Next, we estimate the quantity B from ([9]). For < f < 1, 

p 



B<J2 [ (1 + [di ~ R)+) ln(l + [at - R)+) dx, 

i=l J Bn-l 

while, for 1 < v < 2, to get rid of the nonlinearity given by the factor Qt(a), we rewrite B as 

P r 

/ Qi(a) ln(l + [at - R]+) dx 

4=1 ^n-l 

= J2 / (Qi(o) - + [a - #]+)) ln(l + [ai - rfx 
i=i ■ / B»-i 

+ 22 QiO- + [« - #]+) + [o< - dz. 
i=i •'Sn-i 
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By assumption (j^D the last term is non positive, while (Qi(a) — Qi(l +[a — R}+)) hi(l + [a* — R]+) 
can be estimated via Lemma Finally, from Q and the above, we get 

d ^ f f 

j t J2 (n®(ai-R)dx + 4dJ2 \V x y/l + [ai-R]+\ 2 dx 

1=1 •'^n-l i=l •'Bn 

<C\J2[ {^ + [ai-R]+)H 1 + [ a i-R]+) dx + ^2 f Ha i -R)dx\. 

\i=l JBn-l i= l JB n -l ) 

To conclude, we integrate on (s,t), for — t n < t < 0, and average on s G (— t n -i, — *n)- ■ 

Let us introduce the following easy lemma. 
Lemma 7 Let 

*(z) = (vTTi- i)i { z>0}; Z G K. (10) 
TTie function ^ is nondecreasing, Lipschitz on R, and there exists a constant C > suc/i t/iat 



< Cy^[z), ZER, (11) 
l{*( 2 -M>o} < C2 n iP(z - kn-i), n > 1. (12) 

Proof. To get flUJ), we expand $ and $ both at and +oo to find that there exists C\,Ci > 
with 

V(z) 2 < dinf (z,z 2 ) z>0, 
Hz) > C 2 inf(z,z 2 ) z>0. 

To get (|12|) . first consider z > 2: 
For any n > 1 

2 n tf (z - fc„_0 > - k^) > (1) > (V2 - l)l Wz _ fcn)>0} . 
Next, if z < 2, \l/(z — fc n ) > implies that z > k n and so 2- A; n _i > 2~ n . Hence, for such 2;: 

^(z-kn-!) > [tf (z - - - fc n )] + - fc n ) 

> W(z-k n -.i)-V(z-k n ) 

k n k n _\ 2 



Vl + (^ - fcn-l)+ + V 1 + (* - k n)+ 2^3' 
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We, now, define the sequence U n which plays a key role for the proof of the uniform boundedness 
of the solution. 



U n = sup 2_. / &( a i ~ k n ) dx + 2_. \ \ \V x ^(ai — k n )\ 2 dxds. 

-t„<*<0~ JBn ~{J JQn 

Note that, since \I/ is Lipschitz, 

V x ty(ai - k n ) = V x a/1 + \ai - k n ] + . 
The following nonlinear estimate is a crucial step for establishing our results. 

Lemma 8 There exists a universal constant C ( depending only on A, d, d, P, and N ) such that 

N+2 

Un < CX-1 

for any n > 1. Especially, there exists 5 > (depending only on A, d, d, P, and N) such that if 
Uq < 5 then \im n _ >OQ U n = 0. 

Proof. From Proposition [6] we get the inequality: 

p 

U n < C2 2n J2 l l M 1 + l a i - A; «] + )( 1 + - W dxdT 

p 

■-- 



C2 2n J2[ j l {ai > kn} (l + [ ai -k n } + )^dxdT. 

i—l J J Qn-1 



Hence 

p 



N + 2 



U n <2C2 2n \y2 / lma,-k n )>o}dxdr+\\^(a i -k n )\\ 2 % +2 . (13) 

1 J Qn-1 L 2 ^{Q n -i) 



Next, using ( 1121) . we find that 

p 



i{*(<«-w>fl}^^ 

j = j J J Qn-1 



P 

>N+2 



< (C2 n ) 2 ^ Yl / / - k n-i)?^ dx dr. 

j = X J J Qn-1 



Since ^(a, — fc n ) < \J/(aj — we finally get 

13 



p 

, JV+2 \ — \ 1 1 . 2 iv+2 



U n < 2C2 2n (C2 n ) 2 N ||^( ai - fc n 

f-' ' "7/^(Q„ i) 



Using ( TTTT) . we find that 



By Sobolev imbedding, we find 

\\^(ai - k n _i)\\ 2 2N <cU n -i. 

L 2 (-f„-i,0;L^(B„_i)) 

By interpolation, this gives 

|| ^(cti — /c n _i)|| 2 2(n+2) < cU n -i, with c universal constant. 

(Q n _i) 

Finally, we get 

(iV+2) 

w„ < c n u n _\ 

where C depends only on the constants A as in d, d, P, and N. ■ 

We now discuss the smallness of Uq. 

Lemma 9 For any p > 1 t/iere exists a universal constant C > ( depending only on d, d and p ) 
such that 

p p 

Mo < C(^J ll a illl,P((-3,0)xB(0,3)) + /j W a i\\LP((-3.0)xB(0.3))) ■ 
i=l i=l 

Proof. By definition 

Uq = sup / <&(ai)dx+ / / | V^x/dj + l| 2 da; dr. 

-2<t<0 i=1 Jb{Q,2) J -2 JB(0,2) 

We will use the following facts: 

*(oi) <C(o<(l + |ln(a i )|); (14) 

// jVaj-^aj + l| 2 c?X(ir < / / | V xV /di| 2 da; dr. (15) 
2 JB(0,2) J -2 JB(0,2) 
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The mass conservation yields 

p p p 



d 
dt 



^ / ( (x)ai(t,x) dx = / AACo(ff) a^(t, x) cfe < C ^ 

i=l J B(0,3) i=l JB(Q,3) i=1 



di(t, x) dx : 

■8(0,3) 



where Co > 0, C = 1 on -6(0,2), has bounded second order derivatives and it is supported on 
.6(0,3). Let t G (—2,0). Let r G (— 3,t). We integrate over the time interval (r, £), and then we 
average over r G (—3, —2). Hence, we get 

p . P „ p 

sup / a i (t ) x)dx < CTj / / ai(r,x)dxdr. 

-2<t<0 ~^ JB(Q,2) ~[ J -3 Jb(0,3) 



Similarly, the entropy dissipation yields 

4;Y^ / Co 0*0 Oiln(ai) dx + f (o(x) — i-^^i- rfx 

i=1 </.B(0,3) -/b(0,3) a « 

< / DjACo(^) ajln(aj)dx 

i=1 7^(0,3) 

< C / ai\\n(ai)\ dx. 

i=1 JB(0,3) 

Again we integrate with respect to the time variable. We shall also use the trick 

u\ ln(w)| = u\n(u) — 2uln(u)l < u <i < wln(u) + Cy/u. 

It follows that 

^ ,o 



/ ai\\n(cii)\ dx + d'^^ / / | Vxy/a^ fix dr 

i=l ^(0,2) i=1 J -2 JB(0,2) 

< C I / / ( a i| hi(ai)| + v^i) dxdr\ . 

~i \J-3 Jb(0,3) J 



Combining ffT5]) . (TH]) . and the definition of Uq yields to the desired inequality. 



The proof of Proposition H] follows from Lemma [8] and Lemma [9j 
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4 Duality arguments 



In this section we derive a uniform bound on the solutions a = (a 1; • • -, ap) in a weak norm. We 
will show that this weak norm shrinks through the universal scaling. The proof relies on a nice 
duality argument first used in [26] . 

Let a = (ai, • • •, ap) be a solution to the system §T§ with initial data a = (a°, • • •, a° P ). Then, 
the function p = J2i=i a i * s a solution to the Cauchy problem 



d t p - A (dp) = 0, t > 0, x G 



p 



p°{x) =p(0,x) = J2 a l xe 



(16) 



with 



d[t,x) = 



E i= i a* 

Note that the diffusion coefficient d(t, x) is elliptic. Indeed, it is bounded from above and below: 

< d < d{t, x) < d. (17) 

The equation in f|T6l) seems so to be a nice parabolic equation, except that it is not in the standard 
divergence or non divergence form. In the non divergence form, the equation would provide the 
maximum principle. In the classical divergence form, De Giorgi showed in [7] that such solutions 
are bounded locally in C a . Surprisingly, the behavior of solutions of parabolic equations written 
as (JTBT) is very different. Note that we do not have a priori bounds based on the regularity of d. In 
[26J, Pierre and Schmitt show that any solution of the parabolic equation (1161) . with (1171) (but no 
assumption on the regularity of d), and with regular enough initial values, lies in L 2 ((0,To) x K ) 
for T > 0. However, they also give explicit examples of solutions which blow-up in L p ((0, T ) x W N ) 
for some p > 2. It follows that global regularity of solutions to ([1]) cannot rely only on the equation 
solved by p. 

As usual we denote P(R ) := C(f(lR Ar ). Let us first introduce our weak space. 
Definition 10 We define the space C W (K N ) as the dual space of {p G V(R N ) : Ap G L 1 }: 

C W (R N ) := {/ G V(R N ) : \\f\\ Cw{RN) < oo} 



where 

= sup | < f,if) > | 

!/'GX>(R A, ),||AV|| i i<l 

The following proposition is the main result of this section. 
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Proposition 11 Let a = (ai, • • - ,ap) be a smooth and bounded solution of (QJ) on the time interval 
[0, T] for any T < T (with possible blow-up at T = T ). Then p = Ylf=i a « verifies 

peL»°(0,T ;£ w (^ N )), 

with, 

HpCniU^) < IIpIa^), forallT<T . (18) 
Proof. Let us fix T < T . The proof consists of two steps: 

Step 1: Since a = (a 1; - ■ -,ap) is smooth, the function d is smooth also, except possibly at the 
points (t,x) where p = 0. For small p > 0, let denote a smooth approximation of d verifying 



l^l|L oo ((0,T)xR Ar ) < <^ 



d^(t,x) = d(t,x) if p(t,x)>p. 



Let be a solution to the problem: 



^ - A(d^) = 0, xeR", < t < T, 

p(o,.) = p M (o)ei?(R Ar )nL 1 (R w ). 

We claim that for such a solution, we have for every < t < T: 

To show this claim, consider the dual problem, for any < T < T: 

d t <p + d„A<f) = 0, x e R N , < t < T, 



(19) 



(20) 



(21) 



(22) 



0(T) = # e V(R N ) n L 00 ^). 

The maximum principle, applied to the problem (1221) . gives a uniform bound for the L°° norm of 

0: 

II II £°°(rw) < II^tIU 00 ^)' 

for any < t < T. But, for every < t < T, 



dt 



Pfx{t, x)4>(t, x) dx = 0. 



Then, for any 0^ e Z>(R ) with ||</>f lU 00 ^) — 1> we have 



< 1 and 



p tl (T,x)4> T (x) dx 



p M (0, x)</>(0, x) ofe 



< l|P/ 1 (0)||Li(RiV)||Vi 1 'MI/.' 

< IIpm(°)IIli(k^)- 
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(23) 



niV 



This proves claim (|21[) . 

Step 2: We consider, now, solutions to the following dual problem 

d t (j) + dpAcj) = 0, x e R N , < t < T, 
<p(T) = <p T e V(R N ), 

for any such that 

||A0t||li(RAT) < 1. 

The function 

Pt(t,x) = A0(T-t,x), 0<KT,xGl 
verifies (T20]) with initial value A<pT- Hence, (121]) ensures that 

||A0(t)|Ui(R*) < ||A0 r || L i (RJV) < 1 < t < T. (24) 
But, for any solution a = (oi, • ■ •, ap) of (CQ), p = J^Ci a * ver ifi es (TTOT) and so we have: 

— / p(t,x)4>dx= [ p(t,x)(d — d^)A(f)dx. 

dt J JJiV J R N 

Integrating in time on [0, T], we find 



4>t(x)p(T, x) dx 



< 



0(0, x)p°(x) dx 



/ / p(d — d^)A(j)dx dt 
J Jm 



< IIP°IL 



< \\p°\\c w (RN ) + 2Tp J d. 



+ 2Tpd\\A(P\\ Lao ^ T . L i 



In the second line we have used the definition of the weak norm together with f l24l) . and ( TT9T) . In 
the last line we have used again (1241) . 

Passing to the limit, as p — > 0, for any 0t e V(JSL N ) such that || A^tIIlurjv) < 1, we have 



4>t{x)p(T, x) dx 



< llp u lk 



This implies, thanks to the definition of the weak norm that: 

\\p( T )\\c w (R N ) < l|P°ll£„(R JV )- 

The result follows, since this holds for any T < T . 



As a consequence of this proposition we have the following uniform estimate on the C t 
norm of the total mass. 
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Corollary 12 There exists a universal constant C, such that, for any a = (ax, - ■ -,ap) regular 
and bounded solution to (T7J) on [0,T] x M. N for any T < T (with possible blow-up at T ), we have 
p G L°°(0, T ; C W (R N )), and for any t < T 

\\p(t)\\c w (R N ) < C (||p || Z oo( R iV) + ||p°||li(rJV)) . 

Proof. In the view of Proposition [TTJ we only need to show that the weak norm can be controlled 
by stronger norms. We set 

= r * p° 

where T(x) = N ( N ^_ 2 )uj n \x\*>-z ( u n is the volume of the unit ball in R N ); then, the function solves 

-A0 = p°. 

Since T G L 1 (R N ) + L°°(R N ), there exists a constant C such that 



L oo( R iV) < C(||p°||| L i( M jv) + ||p°||lo 



Therefore, for any if) G V(M N ), with || Aif)\\ L i(uN\ < 1, we have 

p if>\ = | / A(f>if>\ = | / (pAip\ < ||0|| L o O(M iv ) ||A^||Li(RiV) < C(IIp°I IUmk^) + \\p°\\l°° 



We observe that the nonnegativity of the mass p allows us to control stronger norms from this 
weak norm. 

Lemma 13 Let f be a nonnegative function on M. N such that f G L°°(R ) flL 1 (]R Ar ) ; then for any 
compact subset K ofM, N , K CC Mr, there is a positive constant C(K) such that 



\ LHK) <C(K)\\f\\ Cw{RN) . (25) 

Proof. Let K be a compact set in such that K C B(0,r), where B(xq,t) is a ball centered 
at some Xq G K of radius r > . Let G C°°(1R ) be a smooth nonnegative function, such that 
= 1 on K, and = on \j(B(xq, 2r)). Set 

C(tf) = 1^011^). 

Then, the L 1 norm of the Laplacian of f is equal to 1, which implies that 

C(K) 

\\f\\m K) < [ <Kx)f(x)dx<C(K) [ ^f( x )dx<C(K)\\f\\ Cw[RN) . 

■ 

As a consequence of this lemma we have the following corollary. 
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Corollary 14 For any integer P > 0, and any p € (1, j^), there exists a universal constant 
C (depending only on p, P and N) such that the following is true. Given nonnegative functions 
a>i = di(t,x), for 1 < % < P, such that 

p 

p = ^2 a i ^ L°°(—3, 0; jC w (R n )), Vv^G L 2 ((-3,0) x R n ), 1 < i < P, 

i=l 

we have 

a t eL 9 (-3,0; 1/(5(0, 3))), 
for q such that i = 1 — ^ . Moreover, 

_ 2(p-l) 1 

||a»||i«(-3,0;LP(B(0,3))) < C\\ V x v^J £ af(_3 )0 ) XjB W) II Pll L°o(_ 3j0ijCu|(B ,W)) • (26) 

Proof. By Lemma [131 there exists a constant C such that the inequality ([25]) holds 

||p(*,-)||li(b(o,3)) < C||p(t,-)IU ro (K") for - 3<t<0. 
And so, using the nonnegativity of the functions, for any 1 < i < P, we have 

1 1 1 1 (— 3,0;^! (S(0,3))) < C||P| \l°°{-3,0;C w (U. n ))- 

Moreover, by the Sobolev imbedding 

INI l1{ _3,o;l^( R -)) - c \^^ihH(-m^y 

Now, let < 9 < 1, and q defined by 

1 1-9 




N-2 



Then, by standard interpolation, it follows that each a« belongs to the space L q (0, To; L p (B(0, 3))), 
with the bound 

IN|l 9 (- 3 ,0;Lp(B(0,3))) < C(B(0, 3)) || 0* II* (RJV)) II a i II L^(-3,0;LHB(0,3))y ^ 

The inequality f[26l) follows readily by the definition of p. ■ 
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5 Scaling argument 

This section is dedicated to the proof of Theorem [TJ 



The proof is done by a contradiction argument. First, by standard short time existence results, 
we know that there exists a solution for some short time. Next, we assume that the maximal time 
of existence, which we denote by To, is finite. We will show that the L°° norm of a = (ai, • • •, ap) 
is uniformly bounded on (-y , T ) x M N . But this will contradict the fact that the solution blows up 
at T . 



We introduce now the rescaled solutions. Let a = (ai, • • •, ap) be a solution of (OQ). Let 
-y < T < T , x e WL N and < e < a/To/6. Then, we define the rescaled functions 

a £ (s,y) = e^a(e 2 s + T,ey + x ), -— < s < y£R N . (28) 

The function a £ = (a\, ■ ■ ■ , ap) is a solution to a system ([!]) with 



The reaction term Q £ (a) verifies Hypothesis (J2J) , ([3]), (BJ, and ([5]) with the constant A as in ([3]) 
independent on e. The £ w norm of the function p £ rescales by the following identity 



— 2 1 1 I 

|P° IIl°°(— -^-,0;£ U )(R JV )) \\P\\l°°(0,T;Ct, 



This implies, by Corollary [121 that the following bound holds 

\\p £ \\L°°(-3,o;C<»m N )) - C £ ~~ 2 (\\p°\\lhrn) + Hp !!^^)). (29) 
Through the rescaling, from (129]) . we can control the following norm. 

Lemma 15 There exists q > 1 and eq > (depending on d, d, A, P , N , T Q and M ) such that for 
all e < e we have 

p 

||a|||L?((-3,0)xB(0,3)) - 

i=l 

where 5* = S*(q) is the same as in Proposition^ 



Proof. Let < e < a/T /6. We apply Corollary [TD to the rescaled solutions a £ . 
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Each component of G L q (-3, 0; L p (B(0, 3))) for p G (1, ^=5) and g such that I = 1 - -2L; 
Moreover, for each 1 < i < P 

2(P-1) 1 

\\a E i\\L*(-3,0;LP(B(0,3))) < C\\ V^^Of \\^_ 3i q) xR N) IIpIII-^.Q;^^)) ' 

H^x\/°fllL 2 ((-3,0)xR JV ) - Ar |I^V / "illi 2 ((0,To)xR JV )- ( 30 ) 

Hence, together with (J2"5]) . we obtain: 

p 

^ ||l"J(-3,0;Lp(B(0,3))) < C(M), Tq)^^ 
i=l 

where a(p) := ^—{-^1 — 2) + ~(^zt[ — AO and C(M ,T ) is a constant depending on the quantity 
M defined by © and T . 

Since f < 2, the factor (^r — 2) is positive. Thus, there exists p G (1, afz^) such that a := a(p) > 
and 

p 

H a i IU ? ((-3,0)xB(0,3)) < C(M),T )£ Q 

i=l 

with g = (g^j) • This last inequality implies that it is enough to choose £0 < inf 
to conclude the statement for all < e < Eq. 




Next, we apply Proposition @] to the rescaled functions a £ . 



Corollary 16 Let 5* = 5*(q), a and C(Mq,Tq) as in Lemma [731 Then, for all -y < T < Tq, we 

have 
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Proof. Let Xo be any point in W N , -y < T < T . Let a £0 be the rescaled function defined in 
with e = inf |^/^> ( c(a/q t ) ) " j*' Then, by Lemma [TBI we can apply Proposition @] to a e ° to 



obtain: 
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< $>,(x ,T) = e^ 1 J>i°(0,0) < e^P. 

i=l i=l 



Finally, we can now prove Theorem 1. 

Assume by contradiction that the maximal time of existence of a solution a = (oi, • • •, ap), T 
is finite. Then, by Corollary [TBI each a^t, •), for 1 < i < P would be uniformly bounded for all 
-yt < T . In particular, it would be bounded at t = T . Hence, by standard arguments this would 
imply that a £ C°° at t = T . This concludes the proof because it would negate the fact that T is 
the maximal time of existence of a smooth solution a = (ai, • • •, ap). 



A Proof of Theorem 2 



This appendix is dedicated to the proof of the maximum principle in the case P = 2. This is a very 
standard proof. We include it here since we did not find the result in the literature. The system is 
equivalent to the pair of equations 

<9 t ai - diAai = Q( a ), t > 0, x £ R N , 

d t a 2 - d 2 Aa 2 = -Q(a), t > 0, x £ R N , (31) 

a(0,x) = a°(x), xeR N . 

The main remark is that, under Hypothesis (jSJ), we have 

Q(ai, a 2 )(ai — a 2 ) < 0, for any ai, a 2 > 0. (32) 

Let (01,02) be a smooth solution (C 2 ) to ( I3TI) on ([0, T] x R N ), decaying to when \x\ — > 00. For 
any e > 0, we define 

a^(t, x) = ai(t, x) — et, 
a 2 (t, x) = a 2 (t, x) — et. 

Assume that sup(af , a|) attains a local maximum at a point (t e , x e ) £ (0, T] x R^. Let say that its 
value is a\(t £l x £ ). (The proof is similar for the other case.) Especially, this is a local maximum for 
of. Note that 

ft^ — CL 2 — (2i — CL 2 . 

Therefore, from (132p . we have 

<9(a)(t e ,x e ) <0. 
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From the first equation of (T3TT) we have 

(d t a\ - djAafXte, x e ) = Q(a)(t e , x e ) - e < -e < 0. 
This contradicts the fact that a\ attains a local maximum at (t £ ,x £ ). Hence 

sup { a i{t, x), a|(t, x)) < sup (af (t = 0, x), a|(t = 0, x)) = sup a^x)). 

Passing to the limit as e —>■ gives the result. 
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